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This chapter further extends the results obtained in chapters 4 and 5 (from linear 
equation to linear systems). Each algorithm is thoroughly proved and then an example is 
given. 

Five integer number algorithms to solve linear systems are further given. 


Algorithm 1. (Method of substitution) 

(Although simple, this algorithm requires complex computations but is, 
nevertheless, easy to implement into a computer program). 

Some integer number equation are initially solved (which is usually simpler) by 
means of one of the algorithms 4 or 5. (If there is an equation of the system which does 
not have integer systems, then the integer system does not have integer systems, then 
Stop.) The general integer solution of the equation will depend on n — integer number 
parameters (see [5]): 


(P) xX SO sai Le gp 
where all functions f” 
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are linear and have integer number coefficients. 


This general integer number system (p,) is introduced into the other m—1 
equations of the system. We obtain a new system of m—1 equations with n—1 unknown 
variables: 

0, î În, 
which is also to be solved with integer numbers. The procedure is similar. Solving a new 
equation, we obtain its general integer solution: 


WI KS e ), anl, 


where all functions i are linear, with integer number coefficients. (If, along this 


algorithm we come across an equation, which does not have integer solutions, then the 
initial system does not have integer solution. Stop.) 
In the case that all solved equations had integer system at step (j), l<j<r 
(r being of the same rank as the matrix associated to the system) then: 
Î-D 2 FD GY Q) a E n 
(pp) KEP =f (ak = nE, 


f are linear functions and have integer number coefficients. 
J 
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Finally, after r steps, and if all solved equations had integer solutions, we obtain 
the integer solution of one equation with n —r+1 unknown variables. 

The system will have integer solutions if and only if in this last equation has 
integer solutions. 

If it does, let its general integer solution be: 


(r-1) (n EO) () 
P) ko RS kesk 94 = hare li 
where all f® are linear functions with integer number coefficients. 


We’ll present now the reverse procedure as follows. 


r-1) 


We introduce the values of 4; „i. =1,n—r+l1, at step p, in the values of 


ke „i  =l,n-r+2 
from step (p,_,). 
It follows: 


(7-2) _  ¢(r-D (r) (r) (r) (r) (r) (r) 
k; =f J ksk ki” p.k 


| n-r d Sn—-rtl n-r 


(ry g0) (r) 
=8 kesk 


n-r ? 


i_,=ln—r-—1, from which it follows that g” are linear functions with integer 


number coefficients. 

Then follows those (p,_,) in (p,_,) and so on, until we introduce the values 
obtained at step (p,) in those from the step (p,). 

It will follow: 


— ol (pr) (r) 
xi, =i (ki ory Nee 


notation g, k ki =1,n, with all ge. most obviously, linear functions with 
8; 8; y 
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integer number coefficients (the notation was made for simplicity and aesthetical 
aspects). This is, thus, the general integer solution, of the initial system. 


The correctness of algorithm 1. 

The algorithm is finite because it has r steps on the forward way and r -— 1 steps 
on the reverse, (r < +0). Obviously, if one equation of one system does not have (integer 
number) solutions then the system does not have solutions either. 

Writing Sfor the initial system and S, the system resulted from step (p,), 
l<j<r—2, it follows that passing from (p,) to (p;,,) we pass from a system S, to a 


system S,,, equivalent from the point of view of the integer number solution, i.e. 
UD #0 ; Ty? 
k =t,, i; =Lna-j+l, 
which is a particular integer solution of the system S, if and only if 


U) _ p? 
piai 


Tal pi: oe ln-J ? 


is a particular integer solution of the system S,,, where 


SG oa a a= ney 
Hence, their general integer solutions are also equivalent (considering these 
substitutions). Such that, in the end, resolving the initial system S is equivalent with 
solving the equation (of the system consisting of one equation) S_, with integer number 
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coefficients. It follows that the system S has integer number solution if and only if the 
systems S, have integer number solution, 1< j <r—1. 


Example 1. By means of algorithm 1, let us calculate the integer number solution 
of the following system: 


(S) 


Solution: We solve the first integer number equation. We obtain the general 
solution (see [4] or [5]): 


5x — 1y—2z+6w =6 
—4x + 6y-—3z+llw=0 


x=t +2t, 
yt, 

(pi) z=-t,+5t, +3t, -3 
w=t, 


where t, b, t, EZ. 
Substituting in the second, we obtain the system: 


(S,) of = 23t, +21, +9 =0. 

Solving this integer equation we obtain its general integer solution: 
=k, 

(p,) t, =k +2k,+1 


t, =9k +23k, +7 

where k, k, EZ. 

The reverse way. Substituting (p,) in (p,) we obtain: 

x = 3k, +4k, +2 
yak, 
z = 31k, + 79k, +23 
w = 9k, +23k, +7 

where k, k, €Z, which is the general integer solution of the initial system (S). Stop. 


Algorithm 2. 

Input 

A linear system (1) without all a, =0. 

Output 

We decide on the possibility of an integer solution of this system. If it is possible, 
we obtain its general integer solution. 


Method 

1. t=1, h=1, p=1 

2. (A) Divide each equation by the largest co-divisor of the coefficients of the 
unknown variables. If you do not obtain an integer quotient for at least one equation, then 
the system does not have integer solutions. Stop. 


(B) It there is an inequality in the system, then the system does not have integer 
solutions. Stop. 

(C) If repetition of more equations occurs, keep one and if an equation is an 
identity, remove it from the system. 


3. If there is (i, j,) such that 





a, asi then obtain the value of the variable 
x, from equation i,; statement (T,). 


Substitute this statement (where possible) in the other equations of the system and 
in the statement (7_,), (H,) and (P,) for all i, h, and p. Consider t =¢+1, remove 
equation (i, ) from the system. If there is no such a pair, go to step 5. 

4. Does the system (left) have at least one unknown variable? If it does, consider 
the new data and go on to step 2. If it does not, write the general integer solution of the 
system substituting 4, k,,... for all variables from the right term of each expression 
which gives the value of the unknowns of the initial system. Stop. 

5. Calculate 

crt ra; =r mod a, ; 0<|r|<|a,,| 
and determine the indices i, j, j, as well as the r for which this minimum can be 
calculated. (If there are more variables, choose one arbitrarily.) 


a, . . 
6. Write: x, =t, 1— x, , statement (H,). Substitute this statement (where 


ih? 





possible in all the equations of the system and in the statements (7,), (H,) and (P,) for 
all ¢, h, and p. 

7. (A) If az, consider x, =¢,, h:=h+1, and go on to step 2. 

(B) If a=1, then obtain the value of x, from the equation (7); statement 

(P,). Substitute this statement (where possible in the other equations of the 
system and in the relations (7,), (H,) and (P,_,) forall t, h, and p. 

Remove the equation (7) from the system. 

Consider h:=h+1, p= p +1, and go back to step 4. 


The correctness of algorithm 2. Let consider system (1). 


Lemma 1. We consider the algorithm at step 5. Also, let 
M= |r 
Then M#Q@. 
Proof: 
Obviously, M is finite and Mc N“. Then, M has a minimum if and only if 
M +Ø. We suppose, conversely, that M =Ø. Then 
a,, = O(moda, ), Y i Ji jz- 











„a. 


4 =r moda, ; 0<|r|< 


Ba, A A a ie E ARE 


It follows as well that 
a,, = 0(moda, ), V i, ji J- 


That is 











a; =|a,, Vi, joj 
We consider an i, arbitrary but fixed. It is clear that 
(dioan): 4; #9, Vj 


(because the algorithm has passed through the sub-steps 2(B) and 2(C). But, because it 
has also passed through step 3, it follows that 
a, 


Jl, Vj, 

to] 

but as it previously passed through step 2(A), it would result that 
=1, Vj. 


This contradiction shows that the assumption is false. 














Aij 


Lemma 2. Let's consider a,, = r(moda, ). Substitute 


a TT 


a 





X, =t,- 


i Xx 


În 
ig Ja 
in system (A) obtaining system (B). Then 

x, =, jo 


is the particular integer solution of (A) if and only if 
x, =x, j#j, and t, =x}, — 
is the particular integer solution of (B). 


Lemma 3. Let a, z and a, obtained at step 5. 
Then 0 <a, <a, 

Proof: 

It is sufficient to show that a, < 








a,|, V i, 7 because in order to get a,, step 6 is 
obligatory, when the coefficients if the new system are calculated, a, being equal to a 


coefficient form the new system (equality of modules), the coefficient on (i, j,). 








Let a, , with the property 


< 
io Jo Saj; 


a.. 
to Jo 


; there is such an element 





Hence, a, = 








a, | =min a, | . Let a,; with Ca > |a,, 


„j=In 








because la, Fa | is the minimum of the coefficients in the module and not all |a,; 














are equal (conversely, it would result that (a, jor Gay) a, p Vj el,r, the algorithm 
passed through sub-step 2(A) has simplified each equation by the maximal co-divisor of 





its coefficients; hence, it would follow that |a, | =], Vj = Ln , which, again, cannot be 


real because the algorithm also passed through step 3). Out of the coefficients a,, we 


choose one with the property a, #Ma,, there is such an element (contrary, it would 


0 














result (4, jan) la, | but the algorithm has also passed through step 2(A) and it 
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would mean that 





a, | = 1 which contradicts step 3 through which the algorithm has also 
passed). 


Considering qo =[a, a la, [eZ and r=a,, —qa,; EZ, we have 


10] 50 


a <a. We have, thus, obtained an 7, with 


iojo 








=1r,(moda,, ) and 0 salas a,,,|<|a 


oJ. 55 


In| <a, which is in contradiction with the very definition of a,. Thus a, < la, 





, Vij - 


Lemma 4. Algorithm 2 is finite. 

Proof: 

The functioning of the algorithm is meant to transform a linear system of m 
equations and n unknowns into one of m xn, with m <m, n, <n, thus, successively 
into a final linear equation with n—r+1 unknowns (where r is the rank of the 
associated matrix). This equation is solved by means of the same algorithm (which works 
as [5]). The general integer solution of the system will depend on the n—1 integer 
number independent parameters (see [6] — similar properties can be established also the 
general integer solution of the linear system). The reduction of equations occurs at steps 
2, 3 and sub-step 7(B). Step 2 and 3 are obvious and, hence, trivial; they can reduce the 
equation of the system (or even put an end to it) but only under particular conditions. The 
most important case finds its solution at step 7(B), which always reduces one equation of 
the system. As the number of equations is finite we come to solve a single integer number 


equation. We also have to show that the transfer from one system m, xn, to another 
m, is made in a finite interval of time: by steps 5 and 6 permanent substitution of 
variables are made until we to a=1(we to a=1 because, according to lemma 3, all 
a-— s are positive integer numbers and form a strictly decreasing row). 


x Hi 


Theorem of correctness. 

Algorithm 2 correctly calculates the general integer solution of the linear system. 

Proof: 

Algorithm 2 is finite according to lemma 4. Steps 2 and 3 are obvious (see also 
[4], [5]). Their part is to simplify the calculations as much as possible. Step 4 tests the 
finality of the algorithm; the substitution with the parameters k, k,,... has 
systematization and aesthetic reasons. The variables 7, 4, p are counter variables (started 
at step 1) and they are meant to count the statement of the type 7,H,P (numbering 
required by the substitutions at steps 3, 6 and sub-step 7(B); h also counts the new 
(auxiliary) variables introduced in the moment of decomposition of the system. The 
substitution from step 6 does not affect the general integer solution of the system (it 
follows from lemma 2). Lemma 1 shows that at step 5 there is always a, because 
OGzMcN. 

The algorithm performs the transformation of a system m, xn, into another 


m; equivalent to it, preserving the general solution (taking into account, 
however, the substitutions) (see also lemma 2). 


xN. 


i+1? 


Example 2. Calculate the integer solution of: 


12x — 7y+ 9z =12 
— 5y+8z+10w=0 
0z+ Ow=0 
15x + 21z+ 69w =3 
Solution: 
We apply algorithm 2 (we purposely selected an example to be passed through all 
the steps of this algorithm): 
1. t=1, h=1, p=1 
2. (A) The fourth equation becomes 5x + 7z+23w=1 
(B) — 
(C) Equation 3 is removed. 
3. No; go on to step5. 
.a=2 and i=l, j =2, j, =3,and r=2. 
6. z=t +y, the statement (H,). Substituting it in the 
12x —2y +9, =12 
3y+9t,+10w =0 
5x +7y+7t,+23w =1 


Nn 


. a#1 consider z=1¢,,h:=2, and go back to step 2. 


No. Step 5. 
a=1 and i=2, j,=4,j,=2,andr=1. 


Dupa 


. y=t, — 3w, the statement (/7,). Substituting in the system: 
—12x+2t, +9t, —6w =12 
3t, +8t,+ w=0 
5x +7t,+7t,+2w=1 
Substituting it in statement (H,), we obtain: 
z=t,+t, —3w, statement (H,). 
7. w=-3t, — 8t, statement (P). 
Substituting it in the system, we obtain: 
—12x — 20, +57t, =12 
| 5x +h- 9% = 1 
Substituting it in the other statements, we obtain: 
z=10t, +25t, (H,)" 
y=10t, +24t, (H,)" 
h:=3, p:=2, and go back to step 4. 
4. Yes. 
Dies 
3. t, =1—5x+9t,, statement (T). 
Substituting it (where possible) we obtain: 


{1 12x+237t, =—8 (the new system); 
z= 10-50x+115t,, (H )" 
y= 10-50x+1141,, (H,)" 
w=-—3+15x+ 35t, (P)' 
Consider ¢ :=2 go on to step 4. 
4. Yes. Go back to step 2. (From now on algorithm 2 works similarly with that 
from [5], with the only difference that the substitution must also be made in the 
statements obtained up to this point). 
22 
3. No. Go on to step 5. 
5. a=13 (one three) and 7=1, j =2, j, =1, and r=13. 
6. x =t,+2t,, statement (H,). 
After substituting we obtain: 
{-112t, + 134, =—8 (the system) 
z=10-50t,+15t,, (H,)”; 
y =10—50r, +14, (H,)"5 
w=-34+15t, -—5t,, (P)"; 
t= l- 56- t m); 
7. x:=t, h:=4 and go on to step 2. 
Des 
3. No. go on to step 5. 
5. a=5 and i=], j, =1, j, =2 and r=5 
6. t =t, +9t,, statement (H,). 
Substituting it, we obtain : 
St, +13t, =—8 (the system). 
z= 10+85t,+15t,, (HD; 
y= 10+76t,+14t,, (H,)”; 
x= 19t, + 2t, Gia) 
w=-3 -30t — dhs (P)"; 
t= Pa t4 D"; 
7. t =t,;h = 5 and go back to step 2. 
2 = 
3. No. Step 5. 
5. a=2 and i=l, j =2, j, =1 and r=-2. 
6. t, =t; — 3t, statement (H,). After substituting, we obtain: 
St, — 2t, =—8 (the system). 


z=10+85t,-240t,, (H,)"; 
y =10+76t,—214t,, (Hy; 
x= 19f,-55t,, (H,)”; 
w=-3-301,4+85t,, (P)”; 
t,=-1-141,+41t,, (T)"; 


t= 9t,+26t,, (H,)'; 
7. t, =t,,h = 6 and go back to step 2. 
Di 
3. No. Step 5. 


5. a=1 and i=1, j =2, fr. 
6. t, =t; +24, statement (17). After substituting, we obtain: 
t, — 2t, =—8 (the system) 
z =10-395t, —240t,, (H); 
y =10-392t, —2144, (H,)”; 
x= -91t,-55t,, (H3); 
w=-3+140t; +85t, (B); 
t,= 1+68t +41, (TD); 
t= 43 -26t  (H,)"; 
t, = -5t -3t  (H;); 
7. t, =2t,—8 statement (P,). Substituting it in the system we obtain: 0=0. 


Substituting it in the other statements, it follows: 
z = —10304 + 3170 


y = —9184 + 2826 
x =—2371, +728 
w= 3654, —1123 
t, = 1774 —543 

t =112t, +344 
13t + 40 
St, — 16 


statements of no importance. 


E 
t, 


Consider h = 7, p = 3, and go back to step 4. t, EZ 
4. No. The general integer solution of the system is: 
x = —237k, +728 


y =—918k, +2826 
z =1030k, +3170 
w= 365k, -1123 


where k, is an integer number parameter. 
Stop. 


Algorithm 3. 

Input 

A linear system (1) 

Output 

We decide on the possibility of an integer solution of this system. If it is possible, 
we obtain its general integer solution. 


Method 




















1. Solve the system in |”. If it does not have solutions in [|”, it does not have 


solutions in Z” either. Stop. 
2. f=l,t=1, h=l, g=1 
3. Write the value of each main ri X; n the form: 


(E te: ee 


with all qj» de yo Me A, in Z such that all lr, 

















the right term are integer number variables: iie of the secondary variables of the 
system or other new variables introduced with the algorithm). For all 7, we write 
=A: 


E 


4. (E | is dr rX VY, tr =0 where C ; ) are auxiliary integer number 


variables. We remove all the equations (F | ) which are identities. 
5. Does at least one equation (F, :) exists? If it does not, write the general 


integer solution of the system substituting k,,k,,... for all the variables from the right 


term of each expression representing the value of the initial unknowns of the system. 
Stop. 


6. (A) Divide each equation (F ; by the maximal co-divisor of the coefficients 


of their unknowns. If the quotient is not an integer number for at least one i, the system 
does not have integer solutions. Stop. 


(B) Simplify —as in m - all the fractions from the statements (£ ). 

7. Does r „ exist having the absolute value 1? If it does not, go on to step 8. If it 
does, find the value of xX, from the equation (F “ig ) write (T) for this statement, and 
substitute it (where it is possible) in the statements (E ma) (7 >) (H î ) (G,) for all 
i, t, h and g. Remove the equation (F a ) Consider f = f +1, t=t+1, and go back 


to step3. 
8. Calculate 


10 


a = min 
LEI 02 





r 





„7, =r(mod F, ), 0<|r|< 








r.. 
Lp) 


and determine the indices i,, j,, ją as well as the r for which this minimum can be 


obtained. (When there are more variables, choose only one). 


D Q; —r , îi a . 
9. (A) Write x, =z,-—“*—x, , where z, is a new integer variable; statement 


am 
(H, ). 
(B) Substitute the letter (where possible) in the statements 
(£,.) E) (0) (Ea) G,) for alll i t, h and g. 


(C) Consider h:=h+1. 
10. (A) If a #1 go back to step 4. 


(B) If a=1 calculate the value of the variable x, from the equation (F pi B 
relation (G; i) Substitute it (where possible) in the statements (E o (T, ) (H n ). (Ca) 





for all i, t, h, and g. Remove the equation (F,,). Consider g:=g+1, f =f +1 and go 
back to step 3. 


The correctness of algorithm 3 





Lemma 5. Let i be fixed. Then [Ss rX, «| A, (with all 7,, 4, A; n, n, being 
j=m 


integers, n, < n, A, #0 and all xX, being integer variables) can have integer values if 


and only if 7,,,....7,..A, |r. 


in, im i 
Proof: 
The fraction from the lemma can have integer values if and only if there is a 
z €Z such that 
mo, mo, 
e | A =z Dra, —Az+r =0, 
Jan J= 
which is a linear equation. This equation has integer solution © r,,,...,7,,A, |r. 


in, ?°°°? im? i i 


Lemma 6. The algorithm is finite. It is true. The algorithm can stop at steps 1,5 or 
sub-steps 6(A). (It rarely stops at step 1). 

One equation after another are gradually eliminated at step 4 and especially 7 and 
10(B) (F ; ,) - the number of equation is finite. 


If at steps 4 and 7 the elimination of equations may occur only in special cases 
elimination of equations at sub step 10 (B) is always true because, through steps 8 and 
9 we get to a=1 (see [5]) or even lemma 4 (from the correctness of algorithm 2). 
Hence, the algorithm is finite. 


Theorem of Correctness. 
The algorithm 3 correctly calculates the general integer solution of the system (1). 
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Proof: 

The algorithm if finite according to lemma 6. It is obvious that the system does 
not have solution in [| "it does not have in Z” either, because Z” cl)" (step 1). 

The variables f, t, h, g are counter variables and are meant to number the 

















statements of the type E, F, T, Hand G, respectively. They are used to distinguish 
between the statements and make the necessary substitutions (step 2). 

Step 3 is obvious. All coefficients of the unknowns being integers, each main 
variable x, will be written: 


s-[5 cate la 


which can assume the form and conditions required in this step. 
Step 4 is obtained from 3 by writing each fraction equal to an integer variable 
Y,; (this being x, €Z). 


Step 5 is very close to the end. If there is no fraction among all (£ s, ) it means 


that all main variables x, already have values in Z , while the secondary variables of 
the system can be arbitrary in Z , or can be obtained from the statements T, H or G 
(but these have only integer expressions because of their definition and because only 
integer substitutions are made). The second assertion of this step is meant to 
systematize the parameters and renumber; it could be left out but aesthetic reasons 
dictate its presence. According to lemma 5 the step 6(A) is correct. (If a fraction 
depending on certain parameters (integer variables) cannot have values in Z , then the 
main variable which has in the value of its expression such a fraction cannot have 
values in Z either; hence, the system does not have integer system). This 6(A) also 
has a simplifying role. The correctness of step 7, trivial as it is, also results from [4] 
and the steps 8-10 from [5] or even from algorithm 2 (lemma 4). 

Ther initial system is equivalent to the “system” from step 3 (in fact, (£ s, ) as 
well, can be considered a system) Therefore, the general integer solution is preserved 
(the changes of variables do not prejudice it (see [4], [5], and also lemma 2 from the 
correctness of algorithm 2)). From a system m, xn, we form an equivalent system 
m with m,,,<m, and n 


2. 


xn, <n,. This algorithm works similarly to algorithm 


i+] i+] i+] 


Example 3. Employing algorithm 3, find an integer solution of the following 
system: 


3x, +4x, +22x,- 8x, =25 
6x, + 46x, -12x, =2 
4x, +3x,-—x,+9x, = 26 
Solution 
1. Common resolving in 





3 








it follows: 





12 


_ 23x, — 6x; -1 


x, 3 
2 24 
A _Xa+ 2 
llx, +2 
== 
3 








3 Xs 6+% = (E, 2) 
n= -4n e) 


4 X,+2x,-4y, =0 is 
x; —3y,,+2=0 i3 

5. Yes. 

6. — 


7. Yes: |z|=1. Then x; =3y,—2 the statement (7, ). Substituting it in the 
others, we obtain: 





2x,-1 

x, =-7x,+6y,,-4+ 7 (E) 
x,+6y,,-4 

= răi (E) 
3y,,-2+2 

xX, = ~ 12y, +8 + (E) 


Remove the equation (F, 3) 
Consider f = 2, t =2; go back to step 3. 





x, =-7x,+6y,,;-4+ = (E,, :) 
3 NT Ya t5 + = = (E, 2) 
XS —1lyp +8 (E) 
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2X4 oy LO Ei 
4. 
X, +2y3 —4Y =0 Fy» 
5. Yes. 
6. — 


7. Yes |r,|=1. We obtain x, =—2y,,; + 4ypstatement (7, ). Substituting it in 
the others we obtain: 


—4y,,+8y,,—-l ' 
X, =—28 yy +209; oe Ie E, ı 
Ap Yn + Yit5 E, > 
X% = -lly +8 Ezg” 


Remove the equation (F, 2) 
Consider f = 3, t := 3 and go back to step 3. 


3. 
2 2y,,—-1 
x, =—22y,, +30y,, aa BS 
X= Yi + Yon +5 E, 2 
xX = —11y +8 E, 3 
4. 2y,+2y,, +3y,,-1=0 Fey 
5. Yes. 
6. — 
7. No. 
8. a=1 andi, =1, j =31, j, =22, and r=1. 
9. (A) yy =Z —Y (Statement (H,)). 


(B) Substituting it in the others we obtain: 





2y,,+2z,—2y,,—-1 
x, =-22y,, -30z, +30y,,-44+ 28 = AE A PR, 
X= Vat z =y; +5 E 2 ' 


x, =-lly,, +8 É.) 


2y +2z +y, = Ri! 


3,1 


X,=—2y,, +42, -4y,5 T, ' 
(C) Consider h := 2 


10. (B) y =1—2y,; — 2z, statement (G, ). 
Substituting it in the others we obtain: 
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x, =—40y,, —922, +27 E, " 


X= 3y,,+ 3z +4 Eyy" 
a Sally, +8 Ba " 
x, = 6y,,+12z, —4 Ai 
Yo = 2 +32 —l H, ' 


Remove equation (F, L) 
Consider g = 2, f =4 and go back to step 3. 


3. 
x, =—40y,, —922, +27 Ey, 
X= 3y,+ 3z,+4 Ex 
x; == +8 Ea 
4. - 


5. No. The general solution of the initial system is: 
x, = 40k, —92k, +27, from E,, 


X,= 3k, + 3k, +4, from E,, 

x =—I11k, +8, from E,, 
X,= 6k,+12k, —4, from T7," 
w=. 3k —2, from Œ 


where k,,k, EZ. 


Algorithm 4 

Input 

A linear system (1) with not all a, = 0. 

Output 

We decide on the possibility of integer solution of this system. If it is possible, we 
obtain its general integer solution. 


Method 

1. h=1, v=1. 

2. (A) Divide every equation i by the largest co-divisor of the coefficients of the 
unknowns. If the quotient is not an integer for at least one i, then the system 


does not have integer solutions. Stop. 


(B) If there is an inequality in the system, then it does not have integer 
solutions 


(C) In case of repetition, retain only one equation of that kind. 
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(D) Remove all the equations which are identities. 


3. Calculate a = min la, 
UJ 





, Qj 7 O and determine the indices i}, 7, for which this 


minimum can be obtained. (If there are more variables, choose one, at 
random.) 
4. If a#1 goon to step 6. 
If a=1, then: 
(A) Calculate the value of the variable x, from the equation i, note this 


statement (V,). 


(B) Substitute this statement (where possible) in all the equations of the 
system as well as in the statements C) (H k forall v and A. 


(C) Remove the equation i, from the system. 


(D) Consider v:=v+1. 
5. Does at least one equation exist in the system? 
(A) If it does not, write the general integer solution of the system substituting 
k, ko... for all the variables from the right term of each expression 
representing the value of the initial unknowns of the system. 
(B) If it does, considering the new data, go back to step 2. 
6. Write all a, ,, 7 # j) and b, under the form : 





aij = 4; i,j + Mj; > With lr, < la, 





bj g Gain HR > with 


lo 








< 








I, a 
0 


in Jo À 
7. Write x, = = yX; + 9;, +» Statement (H ys 
J#io 
Substitute (where possible) this statement in all the equations of the system as 
well as in the statement ,), (H n) forall v and A. 


8. Consider 
X; 


0 


=t, h=h+]1, 
aij =N J £ Jo 


G0 Ste. ba Sas 


iojo * ~~~ iojo? în ig 


and go back to step 2 


The correctness of Algorithm 4 

This algorithm extends the algorithm from [4] (integer solutions of equations to 
integer solutions of linear systems). The algorithm was thoroughly proved in our previous 
article; the present one introduces a new cycle — having as cycling variable the number of 
equations of system — the rest remaining unchanged, hence, the correctness of algorithm 
4 is obvious. 
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Discussion 

1. The counter variables h and v count the statements H and V , respectively, 
differentiating them (to enable the substitutions); 

2. Step 2 ((A)+(B) + (C)) is trivial and is meant to simplify the calculations (as 
algorithm 2); 

3. Sub-step 5 (A) has aesthetic function (as all the algorithms described). 
Everything else has been proved in the previous chapters (see [4], [5], and 
algorithm 2). 


Example 4. Let us use algorithm 4 to calculate the integer solution of the 
following linear system: 


3% 1 +6x, E2 
4x,+3x, +6x,—-S5x, = 19 
olution 
. h=], v=1 


S 

1 

2. — 

3. a=3andi=l,j=l 
4. 3+1. Goon to step 6. 
6 


. Then, 
-7 =3-(-3)+2 
6=3.2+0 
—2 = 3.0—2 


7. xX, =3x,—-2x,+t, statement (H =). Substituting it in the second equation we 
obtain: 
Ab +3x, +12x, — x, — 5x, =19 
8. x, =t, h=2, ap =0, a, = +2, a,=90, a, = +3, b =-2. 
Go back to step 2. 
2. The equivalent system was written: 
3t, +3x, =-2 
At, +3x, +12x, — x, —5x =19 
3. a=l,i=2,j=4 
4. 1=l 
(A) Then: x, =4t, +3x, +12x, — 5x, —19 statement Y). 
(B) Substituting it in (H, ), we obtain: 
x, =-7t, —6 x, — 21x, +10x +38, (H,) 
(C) Remove the second equation of the system. 
(D) Consider: v:=2. 
5. Yes. Go back to step 2. 


2. The equation +34, +2x, = —2 is left. 
3. a=2 and i=1,j=3 
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4. 2#2, goto step 6. 
6. 
43=42-2-1 
—2 = +2(-1)+0 
7. x, =—2t, +t, —1 statement (H, ). 
Substituting it in (H, J, Y). we obtain: 
x =35t —6 x, -—2lt, +10x, +59 H," 
x, =—20t +3x, +12t, -5x;, -31 V’ 
8. x =t,,h=3, a,=—-1, a; =42, b:=0 , (the others being all = 0). Go back 
to step 2. 
2. The equation —5t, + 27, = 0 was obtained. 
3. a=1,and i=1, j=1 
4. 1=1 
(A) Then ¢, = 27, statement (V, ). 
(B) After substitution, we obtain: 
x =49t, —6x, +10x,+59 H"; 
xXx =-28f6Ć +3x, Dă -31 V"; 
X =—3t, H, '; 
(C) Remove the first equation from the system. 
(D) v=3 
5. No. The general integer solution of the initial system is: 
x, =49k, — 6k, +10k, + 59 


xy = k, 
x, =-3k, -1 
x, =—28k, + 3k, — 5k, — 31 
X; = k, 
where (k, k,,k,)eZ. 
Stop. 
Algorithm 5 
Input 
A linear system (1) 
Output 


We decide on the possibility of an integer solution of this system. If it is 
possible, we obtain its general integer solution. 
Method 
1. We solve the common system in ||”, then it does not have solutions in [l” 
then it does not have solutions in Z” either. Stop. 
2. f=1l,v=1, h=1 
3. Write the value of each main variable x, under the form: 
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x, = GX; =q; f natn | 


j 


with all q;, qi, rj, T, A, from Z such that all HE 





|7; 


(where all x, —S of the right term are integer variables: either from the 





<|AJ, 





A; 7 


secondary variables of the system or the new variables introduced with the 
algorithm). For all 7, we write i= A, 


4. (E, ): dr TX; =T; Ysi +r = 0 where Ca) are auxiliary integer variables. 


Remove al the equations (F ; ) which are identities. 


5. Does at least one equation (F } ,) exist? If it does not, write the general integer 


solution of the system substituting k,, k,,... for all the variables of the right 
number of each expression representing the value of the initial unknowns of 
the system. Stop. 


6. (A) Divide each equation (F, ) by the largest co-divisor of the coefficients of 
their unknowns. If the quotient is an integer for at least one i, then the system 
does not have integer solutions. Stop. 

(B) Simplify — as previously ((A)) all the functions in the relations (E S ). 





#0 , and determine the indices i}, Jọ for which this 


3 Tj 


7. Calculate a = min n]. » 
LJ 


minimum is obtained. 
8. Ifa+1, go onto step 9. 
If a=1, then: 


(A) Calculate the value of the variable x, from the equation (F, ,) write 
V, ) for this statement. 
(B) Substitute this statement (where possible) in the statement (E 7, o) 
(A AR) (H,). for all i, v, and A. 
(C) Remove the equation (£ | P: 
(D) Consider v:=v+1, f= f +1 and go back to step 3. 
9. Write all r ,, j # jo and 7, under the form: 


<|A 











oj =A," dij + oj 2 with n, il> 








Tj = Ain “di +7, , with lr 
10. (A) Write x, = -5 q,,;*; +q, +t, statement H, . 


J#jo 


(B) Substitute this statement (where possible) in all the statements (E 7, o) 


(F,.) C) (H) . 


19 


(C) Consider h := h+ 1 and go back to step 4. 


The correctness of the algorithm is obvious. It consists of the first part of 
algorithm 3 and the end part of algorithm 4. Then, steps 1-6 and their correctness were 
discussed in the case of algorithm 3. The situation is similar with steps 7-10. (After 
calculating the real solution in order to calculate the integer solution, we resorted to the 
procedure from 5 and algorithm 5 was obtained). This means that all these insertions 
were proven previously. 














Example 5 
Using algorithm 5, let us obtain the general integer solution of the system: 
3x FOR, +2x, =0 
4x, — 2x, — 7x, =—l 
Solution 
1. Solving in D ° we obtain: 
—6x, — 2x, 
—2x, +7x,; —l 
p= X +7x; 
4 
2. f=1, v=1, h=1 
—2x, 





3. (Eaa =2x, + 


(E. Jz = e AE = ai. 


4. Fa 22390 


(F ):2x; +3xs — 49. —1=0 
Yes 


i=2 and i, =2, j, =3 
21 
3=2.1+1 
—4=2.-(-2) 
-l1=2-0-1 
10. x, =—x,+2y,, +, statement (H 1): After substitution: 


ne SON 


—2x, 





Eu "x, =2x, —4y,, — 28 + 


3s a =] 


', a 
E2 X B Xs 





F 


1,2 


"tx, + 24, —1=0 
Consider h :=2 and go back to step 4. 
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4. Fa ':-2x,—-3y,, =0 
F ':2t,+x;—-1=0 

5. Yes. 

6. — 


7. a=landi,=2, j,=5 
(A) x; =—2t, +1 statement Y) 
(B) Substituting it, we obtain: 


-2 
A Ae EE A 





Ep "xX, 2h, +1+ yu 
H, ': x, =3t,+1-1+2y,, 

(C) Remove the equation (Fa ). 

(D) Consider v=2, f =2 and go back to step 3. 


=) 
3. Ey, 1 x,=-6t,-4y,+2+—# 


(E,,): x, =-2, +) +1 





4. F, : —2x, eu =0 
5. Yes. 

6. — 

7. a=2 and i,=1,j,=4 
8. 241 

9. —-3=-2-()-1 


10. (A) x, =—y,, +f, statement (H,) 
(B) After substitution, we obtain: 


—2y,, —2t 
E, ': x, = 61, Ay +24 a 
Fa © -yn —2b, =0 
Consider h := 3, and go back to step 4. 
Fy," —Yn—2b =0 
Yes 


eee -® 


a = and i, =1, ją =21 (two, one). 
(A) y; =—2t, statement (V, ). 
(B) After substitution, we obtain: 


(C) Remove the equation F,, . 
(D) Consider v = 3, f = 3 and go back to step 3. 
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3. (E,,): x, =-6t,-4y,, - 2h, +2 


(E; ): xX, =—2t, +N +1 
4. — 
5. No. The general integer solution of the system is: 
x =—ő6k —4k, = 2k; +2, from E;, ; 


x, =—2k, +k, +1, from E,, ; 
X, = 3k, + 2k, =), from H, '; 
44 = 3k, , from A, '; 


x; = —2k; +1, from C€ 


where (k, kk )eZ. 
Stop. 


Note 1. Algorithm 3, 4, and 5 can be applied in the calculation of the integer 
solution of a linear equation. 


Note 2. The algorithms, because of their form, are easily adapted to a computer 
program. 


Note 3. It is up to the reader to decide on which algorithm to use. Good luck! 
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